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A  Stochastic  Model  Providing  a  Rationale 

for  Adjuvant  Chemotherapy* 

James  R.  Thompson  and  Barry  W.  Brown 
Rice  University  and  M.D.  Anderson  Hospital  &  Tumor  Institute 

Abstract  A  model  yielding  the  probability  of  curative  outcome 
for  a  patient  at  the  time  of  tumor  detection  is  presented.  The 
status  of  the  patient  is  determined  by  whether  or  not 
metastases  (distant  spread  of  the  tumor)  have  occurred  and 
whether  any  such  metastases  are  drug  resistant.  If  there  are  no 
metastases,  then  local  excision  is  presumed  curative;  if  there 
are  nonresistant  metastases,  then  local  excision  plus  adjuvant 
chemotherapy  is  presumed  curative;  if  any  metastases  are  drug 
resistant,  there  is  no  cure.  Metastases  and  drug  resistance  arise 
independently  with  intensities  proportional  to  total  tumor  size. 
Over  a  wide  range  of  such  intensities,  the  addition  of  adjuvant 
drug  therapy  yields  a  dramatic  improvement  in  the  probability  of 
cure. 
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Introduction  Goldie,  Coldman  and  Gudauskas  ([4],[5])  have 
explored  the  consequences  of  the  buildup  of  resistance  to  one  or 
more  chemotherapeutic  agents.  They  propose  a  model  in  which 
the  probability  of  the  emergence  of  a  resistant  cell  at  a  given 
time  is  proportional  to  the  size  of  the  primary  tumor  at  that 
instant.  According  to  their  model,  resistance  to  a 
chemotherapeutic  agent  develops  by  mutation,  randomly, 
spontaneously  and  independently  of  any  application  of  that 
agent. 

Another  problem  of  interest  is  that  of  spread  of  the  tumor 
to  other  sites  through  the  process  of  metastatic  progression  [1]. 
Here,  a  cell  may  break  off  from  the  primary  and  move  to  another 
site  in  the  body  setting  up  a  metastatic  clone. 

At  the  time  of  presentation,  a  patient  with  a  solid  tumor  is 
in  one  of  three  fundamental  states: 

(a)  no  metastases 

(b)  metastases,  none  of  which  contain  resistant  cells 

(c)  metastases,  at  least  one  of  which  contains  resistant 

cells. 

Both  simple  excision  of  the  primary  and  the  adjuvant  regime  of 
excision  plus  chemotherapy  will  cure  a  patient  in  state  (a).  A 
patient  in  state  (c)  will  not  be  cured  by  either  regime.  Of 
special  interest  to  us  is  the  probability  that,  at  the  time  of 
discovery  (and  removal)  of  the  primary  tumor  the  patient  is  in 
state  (b).  The  probability  that  a  patient  has  metastases,  all 
nonresistant,  at  the  time  of  presentation,  gives  us  an  indication 
as  to  the  probability  that  a  patient  will  be  cured  by  adjuvant 
therapy  but  not  by  simple  excision  of  the  primary  tumor. 


Discussion  Following  Goldie  and  Coldman  [4],  we  shall  assume 
all  tumors  grow  according  to 

(1)  n(t)  =  e{, 

where  n  is  the  total  number  of  tumor  cells  present. 

The  probability  that  a  nonresistant  tumor  of  size  n  produces  a 
resistant  cell  in  the  time  interval  [t,t+A]  is  given  by 

(2)  Pfresistant  cell  produced  during  [t,t+A] )  =  6  n(t)  a, 
where  1)  represents  the  tendency  to  develop  irreversible 
drug  resistance. 

Similarly,  following  Bartoszyrtski,  Brown  and  Thompson  [1], 
we  shall  assume  that  the  probability  a  tumor  of  size  n(t) 
produces  a  metastasis  in  [t,t+A]  is  given  by 

(3)  Ptmetastasis  in  [t,t+A]  )  =  p  n(t)  a, 

where  p  represents  the  tendency  to  metastasize. 

For  p  and  B  very  small,  we  can  use  n(t)  as  given  simply  by  e*, 
since  here  the  amount  of  tumor  mass  removed  from  the  primary 
to  form  the  metastatic  mass  and/or  the  resistant  mass  is 
negligible.  Furthermore,  we  assume  that  backwards  mutation 
from  resistance  is  negligible. 


Let  us  consider  the  two  events 


(4)  A(N)  =  the  event  that  by  the  time  the  total  tumor  mass 
equals  N  cells  a  nonreslstant  metastasis  develops  in 
which  a  resistant  population  subsequently  develops 
(before  a  total  tumor  mass  of  N). 

(5)  B(N)  =  the  event  that  by  the  time  the  total  tumor  mass 
equals  N  cells  a  resistant  clone  develops  from  which  a 
metastasis  develops  (before  a  total  tumor  mass  of  N). 

We  shall  seek  to  compute  P(AC(N))  and  P(BC(N)).  Now, 


(6)  P(metastasls  occurs  in  [t,t+A]  followed  by  a  resistant 

subclone  before  T)  =|in[1-exp{-j0T  tJ3eTdx}]A 

=  |iel[1  -  eBexp(-BN/n)], 
where  T  =  In  (N). 


But  then, 

(7)  P(AC(N))  =  expl-l^iie^l-e^xpf-BN/nMdt] 

=  exp[|i-|ieBe“BN  +iiBeBN(Ei(-0)  -Ei(-BN)}] 
where  the  exponential  integral  Ei(x)  is  defined  for 
negative  x  by 

(8)  Ei(x)  =  -  J  -  e“l/t  dt. 

'-x 

Similarly,  we  obtain 
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(9)  P[BC(N)]  =  explB-Be^e^+pBe^N  (Ei(-p)-Ei(-uN)}  ] . 


Thus,  we  can  compute  the  probability  that  no  resistant 
metastases  have  been  formed  by  a  total  tumor  size  of  N  via 

(10)  P(no  resistant  metastases)  =  P(AC(N))  P(BC(N)). 

We  have  been  able  to  obtain  (7)  and  (9),  exact  solutions  to  the 
special  two  stage  branching  processes  defined  by  hypotheses 
(l)-(3),  (generally  deemed  nontractable  and  requiring 
approximation)  by  exploiting  the  special  structure  of  the  case  of 
interest  to  us  here:  namely,  the  nonappearance  of  any  second 
stage  events.  The  exact  expression  in  (10)  is  easily  computed 
using,  for  example,  the  IMSL  routine  MliDEI  [6]  for  computing 
exponential  integrals.  However  for  the  magnitudes  of  N  (very 
large)  and  p  and  B  (very  small)  in  the  present  application,  we 
obtain  essentially  the  same  result  (using  approximations  [7]  for 
small  arguments  of  Ei(x)  )  with 


(11)  P(no  resistant  metastases  by  N)  =(pB82)wBN 
where  6  =  1.781072  is  e  raised  to  Euler's  constant. 

The  probability  of  no  metastases  at  all  by  tumor  size  N  is  given 
simply  by 

(12)  P(no  metastases  by  N)  =  exp[-j0Tpeldt] 

=  exp[-p(N-1)]. 
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Thus,  the  probability  a  patient  with  a  tumor  of  size  N  is  curable 
by  adjuvant  therapy  but  not  by  simple  excision  of  the  primary 
tumor  alone  is  given  by 

(13)  P(metastases,  none  of  them  resistant  by  N)  = 

(pBa2)1*6"  -  exp[-p(N-1)J . 

In  Figure  1,  we  show  the  probability  of  nonoccurrence  of 
metastases  versus  log  (pN).  In  Figure  2,  we  show  the 
probability  of  nonoccurrence  of  resistant  metastases  versus 
log(N)  for  various  values  of  log(pB).  Typical  values  for  J3  for 
mammalian  cells  are  10~6  to  I0~4  [4].  A  typical  tumor  size  at 
time  of  detection  is  1010  cells  (roughly  10  cubic  centimeters). 
In  Figure  3,  we  show  the  probability  there  exist  metastases  but 
no  resistant  metastases  at  a  tumor  mass  of  1010  for  various  B 
and  p  values.  We  note,  for  example,  that  for  B=10~6-5  to  10~10, 
the  probability  a  patient  presents  with  a  condition  curable  by 
adjuvant  therapy  but  not  by  simple  excision  of  the  primary  is  at 
least  40%.  Similar  results  hold  for  early  detection  ( 1 09  cells) 
as  shown  in  Figure  4.  This  wide  range  of  p  values  includes  that 
reported  for  breast  cancer  by  Bartoszyriski,  Brown  and  Thompson 

m. 

An  alternative  procedure  to  that  above  can  be  based  on  an 
approximation  used  by  Goldie,  Coldman  and  Gudauskas  in  another 
context  [5].  Using  the  axioms  (1)-(3),  letting  R(t)  be  the  number 
of  resistant  cells,  we  may  write 

( 1 4)  dE(R(t))/dt  =  E(R)  +  Bet  =  E(R)  *  Bn(t) 
with  solution 
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(15)  F(R)  =  Bn  In  n 

Then,  the  probability  for  no  cells  which  develop  first  resistance 
and  then  metastasize  is  approximated  by 

(16)  Pr[no  metastases  thrown  off  by  resistant  clones!  « 

exp(-J0Tu  E(R)  dx)=  exp[-pB(  1  +N{lnN- 1 ))]. 

Similarly,  we  have  that 

(17)  Pr[no  metastasis  which  develops  resistance]  ~ 
exp{-|iB(  1  +NClnN- 1 })]. 

Then  the  G-C-G  approximation  to  (1 1)  is  given  by 

(18)  Prfno  resistant  metastasis  by  N]  ~ 
expl-2iiB(1+N{1nN-1})]. 

For  probabilities  greater  than  .5,  the  G-C-G  approximation  is 
very  close  to  (1 1). 

It  has  been  reported  [8]  that  in  vivo  B  values  may  actually 
be  a  hundred  times  greater  than  in  vitro  values.  For  B  values  of 
.01  or  greater,  the  assumption  that  the  nonresistant  population 
at  time  t  is  given  by  e*  is  inappropriate.  Accordingly,  we 
consider  below  an  approximation  based  on  a  modification  of  the 
G-C-G  argument. 

(19)  d(E(R))/dt  *  E(R)  +  B[  el  -R] 
has  the  solution 

(20)  E(R)  *  e*  -  et(1_B). 

Then, 

(21)  Prfno  metastasis  thrown  off  by  resistant  clone  by  NJ  ~ 
expl-u(N-l)  +p/(1-B)  {N1_B-1)1. 

Similarly, 

(22)  Pr[no  metastasis  thrown  off  which  develops  resistant 
clone  by  N]  *  exp[-B(N-1 )  +B/(l-p)  (N1-^-!)]. 


Summarizing, 

(23)  Prfno  resistant  metastasis  by  total  tumor  mass  N]  * 
exp[-(B+p)(N-l)  +  p/(l-B)(N(,_6)-1)  +B/(  1  -  p)(N  (,~  *  )  -1)1. 

By  differentiating  (21)  and  (22)  with  respect  to  N,  we  note  that 
for  B  very  large  relative  to  p,  the  chances  are  that  any 
formation  of  a  resistant  metastasis  after  tumor  discovery 
would  most  likely  be  the  result  of  spread  from  a  resistant  clone 
in  the  primary.  Here,  the  standard  protocol  of  removing  the 
primary  before  beginning  chemotherapy  would  be  indicated.  On 
the  other  hand,  if  the  force  of  metastasis  if  much  stronger  than 
that  of  mutation  to  resistance,  then  chemotherapy  might 
appropriately  precede  surgical  intervention. 

In  the  case  of  breast  cancer  which  has  metastasized  at 
least  to  local  nodes,  it  has  been  reported  [2]  that  the  use  of 
adjuvant  chemotherapy  decreases  disease  mortality  by  as  much 
as  54%  when  compared  to  surgery  alone.  In  order  to  estimate  p 
and  B  clinically,  we  need  randomized  trials  on  tumors  (which 
have  exhibited  no  metastases  at  presentation)  using  surgery 
followed  by  adjuvant  chemotherapy  compared  with  surgery  alone. 
Such  a  clinical  data  base  is  not  currently  available.  However, 
since  we  are,  at  this  stage,  really  seeking  rough  estimates  of  p 
and  B,  animal  experiments  may  be  appropriate. 

We  examine  below  how  such  experiments  might  be  used  to 
estimate  p  and  B.  Let  us  suppose  we  have  stratified  our  data  by 
tumor  size  at  presentation.  Consider  the  lO10  primary  cell 

stratum.  Suppose  that  the  control  group  (surgical  excision  only) 
exHitms  a  cute  rate  at  9»  ana  that  the  adjuvant  therapy  group 

exhibits  a  95%  cure  rate.  Then  we  can  estimate  p  from  (12): 

(24)  expl-p(N-l)]  =  expl-  p  (  10to-1)]  =  .05. 
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This  yields  a  ji  value  of  3  *10  10 
Next,  we  can  estimate  B  from  (11): 

(25)  (Aju2)^  =(  3  *10~,0£  *2)  3  *  =  95. 

This  gives  (Using  Newton's  Method)  £  =  .0006.  (The  same 
estimate  is  obtained  by  the  use  of  (23).  ) 

Although  we  are  here  essentially  concerned  only  with 
obtaining  rough  estimates  of  n  and  £,  it  is  clear  that  a  number 
of  resampling  techniques  (e.g.,  the  jackknife  or  the  bootstrap 
1 31),  can  be  used  to  determine  the  variability  of  the  estimates. 
Let  us  suppose,  for  example,  that  we  have  Nj  individuals  in  the 
control  group  of  whom  n}  are  cured  and  N2  individuals  in  the 
adjuvant  group  of  whom  n2  are  cured.  Using  the  coding  of  1  for 
a  cure  and  0  for  a  noncure,  we  repeatedly  sample  (say  M  times, 
where  M  is  very  large)  N,  individuals  with  probability  of 
success  nt/Nt  and  N2  individuals  with  probability  of  success 

n2/N2.  For  each  such  j'th  sampling  we  obtain  Aj  and  as 
above.  Then  we  have  as  ready  bootstrap  estimates  for  Var(ji), 
Var(#),and  Cov(A,i),  2(  jij-ji)2/M,  2({j-#)2/M  and 

S(ji|-ji)(,pj-£)/M,  respectively. 

Conclusions  We  have  developed  a  model  which  gives  a 
framework  in  which  the  potential  benefits  of  adjuvant 
chemotherapy  can  be  considered.  There  seems  to  be  a  strong 
indication  from  the  model  that  chemotherapy,  in  conjunction 
with  surgery,  is  frequently  preferred  to  surgery  alone,  not 
simply  on  considerations  of  palliation,  but  on  considerations  of 


cure. 
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